Domain wall switching: optimizing the energy 

landscape 

Zhihong Lu 

MINT Center and Department of Physics and Astronomy 
University of Alabama, Tuscaloosa AL 35487-0209 
Email: zlu@mint.ua.edu 
P. B. Visscher 
MINT Center and Department of Physics and Astronomy 
University of Alabama, Tuscaloosa AL 35487-0209 
Email: visscher@ua.edu 
W. H. Butler 

MINT Center and Department of Physics and Astronomy 
University of Alabama, Tuscaloosa AL 35487-0209 
Email: wbutler@mint.ua.edu 



Abstract — It has recently been suggested that exchange spring 
media offer a way to increase media density without causing 
thermal instability (superparamagnetism), by using a hard and 
a soft layer coupled by exchange. Victora has suggested a figure 
of merit ^ = 2Et/ ^DmsHsw, the ratio of the energy barrier 
to that of a Stoner-Wohlfarth system with the same switching 
field, which is 1 for a Stoner-Wohlfarth (coherently switching) 
particle and 2 for an optimal two-layer composite medium. 
A number of theoretical approaches have been used for this 
problem (e.g., various numbers of coupled Stoner-Wohlfarth 
layers and continuum micromagnetics). In this paper we show 
that many of these approaches can be regarded as special cases 
or approximations to a variational formulation of the problem, 
in which the energy is minimized for fixed magnetization. The 
results can be easily visualized in terms of a plot of the energy 
i5 as a function of magnetic moment m^, in which both the 
switching field [the maximum slope of i5(mz)] and the stability 
(determined by the energy barrier Ai5) are geometrically visible. 
In this formulation we can prove a rigorous limit on the figure 
of merit ^, which can be no higher than 4. We also show that a 
quadratic anistropy suggested by Suess et al comes very close to 
this limit. 

I. Introduction 

Recently the concept of an exchange-spring medium[l], 
[2], [3] whose grains have a soft and a hard layer has 
been generalized [4] to a system with a continuously- varying 
anisotropy. Various models of this system have been explored 
- the purpose of this paper is to show that the relationships 
between these can be easily visualized by using a variational 
formulation of the problem. 

In an exchange-spring medium, we want to minimize the 
switching field, but of course we can make this as small as 
we want by using a very small anisotropy, and the medium 
will be superparamagnetic (thermally unstable) and useless. To 
make comparisons between media, we must hold something 
constant to maintain stability. Often what is held constant is the 
anisotropy field (the coercivity IKjM^ of the hardest layer. 
However, this hardest layer might be very thin and have little 



effect on the overall coercivity. In a practical application the 
quantity it is most important to hold constant is the overall 
energy barrier to switching, which we will assume determines 
the thermal stability of the medium. Victora[l] introduced a 
figure of merit for this purpose, ^ = 2Eb/ ^qitIsHsw (here Eb 
and TTis are the barrier energy and saturation magnetic moment 
per unit area) for which we will prove a rigorous bound in Sec. 

uni 

To do this, we develop a variational formulation in which we 
describe the switching behavior in terms of a function E{raz), 
the energy per unit area as a function of the magnetic moment 
per unit area. This turns out to be a very useful way of thinking 
about switching problems. 

II. Model 

We consider a one-dimensional model, in which the magne- 
tization M(z) is a function only of one variable z (independent 
of X and y). We will allow the anisotropy K{z), exchange 
constant A{z), and saturation magnetization Ms{z) to vary ar- 
bitrarily with z. Since we will do computations with a discrete 
approximation to this continuum model (which approaches 
the continuum model as the cell size 0), we will write 
the energy in a discrete form. It has cells labeled by i, with 
magnetization vectors M^. In the quasistatic energy minima 
we will consider, these vectors will lie in a plane, so they 
can be described by giving the angle 6i of the magnetization 
relative to the long axis of the grain (the z axis): 

The energy (per unit area in the xy plane) E of our system 
is then given in terms of the values of K and M at each cell 
(and A between each neighboring pair of cells) by 

E = ^ aiKi s\T? 9i + ^ ai — ^'^^ cos{9i+i — cos Oi) + 

N 

'^QiHoMiH COS 9i 

i=l 



where Ki is the (perpendicular) anisotropy at the center of 
cell i, Qi is the length of cell i, i+i is the distance between 
cells i and i + i+i is the continuum exchange parameter 
evaluated between these cells [5], and H is the external field 
(assumed along z). For simplicity, we do not consider mag- 
netostatic energy here - in similar systems, micromagnetic 
simulation has shown that this affects the coercivity by only 
a few percent. 

We consider here quasistatic switching - we assume that 
we vary H in such a way that the system is always at a 
relative minimum (with respect to the M^s) of the energy. 
More precisely, we assume H is very slightly above this value, 
so that the magnetic moment 
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(2) 



increases slowly, and we consider the limit in which the rate 
of increase approaches zero. Note that this is never true in a 
real switching event - after a domain wall has traversed most 
of the sample, it would require reversing H to keep the system 
quasistatic. However, by this time it is irrelevant whether the 
system remains quasistatic (it will finish switching in either 
case) and in the initial stages the quasistatic assumption is 
often reasonable. 

It would appear that to find the quasistatic switching trajec- 
tory, in which H varies with time, we would need to minimize 
a function E{9i, 62, - --dN, H) of a large number of variables 
61,62, ■■■On, for each value of H independently. However, 
there is a way around this. We can choose some coordinate 
in the space of diS (we choose the longitudinal component of 
the magnetic moment, m^, for reasons apparent below) and 
first minimize £{61,62, ■■■On, H) for fixed m^, obtaining a 
function (the constrained minimum energy) E{mz,H)^ Then 
we can minimize E{mz,H) with respect to m^, obtaining 
the same relative minimum E{H) we would have obtained 
by unconstrained minimization. [Note that there may be more 
than one relative minimum, so we should call this Ej{H) 
where j indexes the minima, but we will omit this index 
for simplicity.] The advantage of this apparently-circuitous 
method of finding the minimum is that the configuration 
minimizing E{6i,02, ■■■On ,H) is actually independent of H\ 
This is apparent from Eq. ([T]i above, since the only dependence 
on H is the Zeeman term fiQUizH, which is a constant when 
rriz is held fixed. The result is that we need only compute the 
constrained minimum energy at iJ = 0, and it is given at any 
other field H by 



E{mz,H) = £'(777,2, 0) - IJ-om^H 



(3) 



Furthermore, this energy is minimized at a particular H by 
setting dE{mz, H) /dniz = 0, so the field necessary to hold 
rriz constant is given by 

dE{m„0) 



fioH 



dm. 



(4) 



We conclude that everything we need to know about the 
system (the coercivity and energy barrier) is contained in the 



function E{mz), the minimum energy at fixed magnetic mo- 
ment mz and zero field. This result is very general. Although 
we motivated it above by considering domain-wall switching, 
it describes Stoner-Wohlfarth (S-W) switching as well. This 
is the limit in which K, A, and Ms are uniform and A is 
large so M(z) is uniform. The S-W energy (per unit area, of a 
grain of length L) is just E = KL sin^ = KL{1 — ml /ml) 
(here the saturation moment per unit area is ms = MgL) so 
the E{mz) plot is a parabola, as shown in Fig. [T] Note that 




mz/ms 

Fig. 1. Energy landscape E{mz) for a Stoner-Wohlfarth particle. 

the slope (the field necessary to switch, Eq. |4|i is exactly the 
Stoner-Wohlfarth switching field Hsw — 2K/Ms, as expected. 
[In general, the coercivity is the maximum value of the slope.] 

Another virtue of the function E{mz) is that it has exactly 
the same interpretation for a particle with a lower exchange 
constant A as for a S-W (high-^) particle. The behavior 
depends only on the dimensionless parameter x = A/KL'^, 
which is the square of the ratio of the exchange length to the 
particle length L. With low x, switching takes place through 
domain wall motion, and the energy barrier is approximately 
the domain wall energy. This can be calculated analytically 
in an infinite system (we will refer to this as the thin- 
wall approximation, because it is valid when the wall is 
far from the system boundary and the material properties 
vary only slightly through the wall) - the thin-wall energy 
is 4(v4ii')2 = 4iirLx2. We have developed a numerical 
minimization program for computing E{mz) for an arbitrary 
K{z), and the result for a uniform K{z) is shown in Fig. |2] 
It can be seen that the energy is indeed constant when mz 
is far from its limiting values ±ms = ±MsL, and equal to 
A{AK)2. The slope of E{mz) at the ends is just the domain 
wall nucleation field. The behavior of the magnetization profile 
at various times during switching is shown in Fig. [3] 

Suess et al[6] have noted that in the thin-wall approxima- 
tion, the pinning field should remain constant if we choose 



Kiz) 



. In our E{mz) formulation, this means the slope 



should be nearly constant. This turns out to be remarkably 
nearly true numerically, except near the hard end, as shown in 

Fig. a 
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Fig. 2. Energy landscape E(m2) for a particle with small exchange 
parameter x = 0.0076. 
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Fig. 3. Magnetization profiles (angle vs. position z) for the particle whose 
energy landscape is shown in Fig. |2] Labels (a), (b), ... correspond to specific 
values of niz shown in Fig. |2] 



III. Rigorous bound on coercivity figure of merit 

Our E{mz) formulation allows us to prove a completely 
general (within the assumptions: ID, quasistatic) result, which 
is clear geometrically from the E{mz) graph. If we fix the 
vertical height (the zero-field barrier E},) and the horizontal 
extent {2ms) the minimum possible coercivity (coercivity = 
maximum slope) is obtained by a straight line, whose slope 
must be pioH = Eb/2ms. In terms of the figure of merit, this 
means ^ < 4. 

Another way of stating this result is that the coercivity of 
any graded medium cannot be less than 1/4 of the coercivity of 
a Stoner-Wohlfarth particle (assuming that the latter switches 
coherently) of the same magnetic moment and energy barrier 
[If Ms is constant, fixing the moment is the same as fixing 
the length L.] Note that the K{z) oc case (Fig. |4]l gives 
^ = 3.23, close to the theoretical limit, which it approaches 
as X — > 0. 

Note that in this paper we consider only fields along the 
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Fig. 4. Energy landscape E{mz) for the case K(z) = {Ki^/ L'^)z^ , which 
would be exactly linear in the thin-wall approximation. We have used Kh = 
6 X lO'^J/m^, L = 22.2 nm, A = 1.0 X 10"" J/m, so that the dimensionless 
exchange parameter AjKyi^l? = 0.00338. 



easy axis. Obviously it is worth considering how transverse 
fields might be useful in switching, since it is known that 
by using a field at 45° from the axis the Stoner-Wohlfarth 
switching field is decreased by a factor of 2, so the figure of 
merit ^ becomes 2. Also, it is likely that the nucleation of 
a domain wall (which initially requires transverse twisting of 
the magnetization) can be assisted by a transverse field, so 
the figure of merit might increase slightly above 4. We should 
also note that the limit we have established assumes fixed grain 
length - it may be possible to decrease switching fields beyond 
the factor of 4 because graded media may make it possible to 
use longer grains without encountering complicated switching 
modes such as vortices. 

IV. Conclusion 

We have shown a general bound on the figure of merit of 
a graded-anisotropy medium, and that this bound (4) is very 
nearly achieved by K{z) oc z^. Because the usefulness of 
such a medium depends on its thermal stability as well as its 
coercivity, and because of the complex switching mechanism 
the zero-field switching rate is not completely determined 
by the energy barrier, an important remaining problem is 
the more precise calculation of this rate. Although brute 
force micromagnetic simulation of such slow switching is 
not practical, work is under way on accelerated sampling 
techniques for solving this problem[7]. 
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